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Abstract 



We review the status of the electric dipole moment (EDM) of neutron in the Standard 
Model (SM). The contributions of the strong and electroweak interactions are discussed 
separately. In each case the structure of the Lagrangian and the sources of CP violation 
are specified, and calculational details are given subsequently. These two contributions 
to the neutron EDM exist in any extension of the SM including supersymmctry, two- 
doublet models as well as models with more than three generations of fermions. We do 
not discuss the EDM in such extensions; however, we briefly summarize their predictions 
with a detailed account of the related literature. 
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Introduction 



The electric dipole moment (EDM), d, of a classical charge distribution p{x) is given by 

<r=/d-'rfp(f). (1) 

which is a polar vector. It is known that the elementary particles have no intrinsic vector 
quantity other than their spin. Therefore, their electric (d) as well as magnetic dipole 
(/I) moments are expected to be proportional to their spins (S). So under discrete space- 
time transformations, these moments and the spin must have the same transformation 
properties. 

Despite their similarity concerning the particle spin, the electric and magnetic mo- 
ments have quite different characteristics. This can be understood by studying the par- 
ticle under electromagnetic field. Depicted in Fig. 1 is a particle (grey blob) with spin S 
under electric (E) and magnetic (B) fields. Consider the upper part of the figure. Ini- 
tially, particle is in a frame where B\\S. Under under time-reversal (T) operation (where 
B ~ B and S ^ — S) they still remain parallel. Therefore, the corresponding mo- 
ment of interaction, that is, the magnetic moment (/I) implies no violation of the discrete 
symmetry T= CP, assuming CPT is a good symmetry of the Nature. 

Consider now the lower part of the same figure where initial configuration is similar 
to that of the upper part; -^115*. However, under a T operation E ^ E and S — S, 
and thus, the corresponding interaction term in the Lagrangian violates the CP symme- 
try. Therefore, unlike the magnetic moment discussed above, the electric dipole moment 
implies the violation of CP invariance in Nature. It has long been known that CP is not 
a conserved symmetry of the Nature. Indeed, the CP-odd long-living neutral kaon is 
known to disassociate into the CP-even final state consisting of two neutral pions tt'^tt" 
[|I], 0. Therefore, CP is not a respected symmetry at all so all particles with spin, must 
have an EDM at some level. 

At the field theory level, the P and T violating interaction between a spin 1/2 fermion 
field ipi^) and the photon can be written as: 

- D{e)i^ (g + k/2) (T^.75A;> (g - k/2) A''{k) (2) 

where k is the momentum taken away by the photon. One notes that when 75 is removed 
from this operator the corresponding quantity becomes the magnetic moment interaction 
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Figure 1: Illustration of S ■ B (up) and S ■ E (down) interactions under time-reversal (T 
= CP). The former refers to the magnetic dipole moment of the particle whereas latter to 
the electric dipole moment. Clearly a finite EDM implies the existence of the CP-violating 
operators in the underlying theory. 



which does not imply the discrete symmetries. The value of the form factor D{k'^) for 
on-shell photon (/c^ = 0) is the sought expression for the EDM of the fermion ip: 

d^ = -D{k^ = 0) (3) 

Of all the baryons the EDM of the neutron Q dn has been of interest to the physicists for 
a long time. In 1950, Purcell and Ramsey first considered the problem of existence 
of the neutron EDM. However, at that time, it was assumed that the physical world is 
invariant under P inversion, and in order to have a nonzero dn, Purcell and Ramsey had 
to construct a rather unconventional P even EDM of neutron. They gave the first limit 
on the EDM of the neutron: 

dn=\dn\< 10~^3-3.10~^^ e - cm . (4) 

In 1957, Landau observed that a non- vanishing EDM of the neutron was a signal of P 
and T violation. About that time it was discovered that P invariance is not a strict law 
of nature |0 as was experimentally confirmed [0. EDM of neutron violates both parity 
and time reversal invariance. 

Even after the parity was found to be violated by the weak interactions it was be- 
lieved that EDM of the neutron was time-reversal invariant and by virtue of CPT theorem 
was also invariant under CP symmetry. In 1964, when the CP violation was observed in 

^The neutron is most convenient object for precision measurements of the EDMs, owing to vanishing 
of its electric charge, its sufficient stabihty, and the possibihty of using beams of uhracold particles. 
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kaon system the subject of EDM of neutron has become of particular importance for 
both theoretical and experimental interests. 

Various methods have been invoked for the experimental investigation of the EDMs 
of elementary particles (for details see, for example, 0). So far only the upper limits on 
the EDM of different particles have been obtained [§. 

The detection of a non-zero EDM for neutron would have extremely important im- 
plications. It would be the first observation of CP violation outside the kaon system. So, 
it could give us a new understanding of the CP violation mechanism and of the early uni- 
verse cosmology where it is generally believed that CP violation gave rise to the observed 
baryon-antibaryon asymmetry of the universe |10 . 



Since the experimental discovery of CP noninvariance in Nature, many diverse mech- 
anisms of the CP violation have been proposed and several estimates of the EDM of 
neutron have been obtained. A comprehensive review of the EDM of neutron in various 



models of CP violation can be found in 11 



In this dissertation we shall consider the CP violation and hence the neutron EDM 
in the frame work of SU{3)c x SU{2)l x f/(l)y theory of electroweak and strong interac- 
tions, i.e., The Standard Model. In the standard model there are two sources of CP 
violation namely: 

1. The first source of CP violation in the standard model is related to the properties 
of strong interaction theory described by Quantum Chromodynamics (QCD). It has long 
been realized that, the nontrivial structure of vacuum in QCD and the existence of in- 
stantons [0 in the non-Abelian gauge theory make it necessary to add to the Lagrangian 



the so-called 6'-term, which violates P and CP symmetries. This term for 9 ~ calO{l) 
gives an exceedingly large value for the neutron EDM which can be made consistent with 
the experimental bounds [|] iiO < 10~^. Several methods have been employed to evaluate 
the effect of CP violation on neutron EDM. The first investigation on this subject was 



carried out by Baluni [Q. In Chapter I we will discuss sources of CP violation in the 



strong interactions together with a calculation of the neutron EDM. 

2. The second source of CP violation in the SM appears in the electroweak sector where 
the mixing matrix of quarks VkM) appearing in the charged-current interactions, possess 
a complex phase which generates CP-violating observables. For three generations of the 
chiral fermions this matrix can be parametrized by three Euler angles and a single phase. 



5cKM |T6[. For two generations, for instance, there is no CP violation because the 
quark mixing matrix can always be made real. The analysis of the EDM of the neutron in 
the electroweak sector is presented in Chapter II, where we have adopted the calculational 
techniques given by Nanopoulous et al 
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Chapter 1 

The Neutron EDM: Strong 
Interactions 



In this chapter we will study first the sources of CP violation in strong interactions. For 
this purpose we will work out the nonperturbative and topological properties of the QCD 
vacuum together with its symmetries. Progressively, we will construct the representation- 
independent CP violation Lagrangian which is the mere source of CP-violating quantities 
in the hadron spectrum. Finally we will compute the EDM of neutron and compare the 
theoretical estimate with experimental bounds. 



1.1 Symmetries of the QCD Lagrangian 



The quantum chromo dynamics (QCD) is the gauge theory of quarks and gluons, which 
are believed to be constituents of the hadrons. The total action density of the theory can 
be divided into two distinct parts, 



^QCD — ^0 + 

where Cq describes the quarks (color triplets: ip"" , 
, a = 1, ■ ■ ■ , 8) together with their interactions 



(1.1) 

1, 2, 3) and gluons (color octets: 



^9. 



Tr 



^R^g^L + h.C 



;i-2) 



where gs is the gauge (here color SU{3)) coupling constant. In this expression the gluon 
field strength tensor reads as 



d^A^ — di,A^ + i\A^^ Ap 



where A^ = y4°A", and 3x3 anti-hermitian matrices A" (a = 1, 
generators satisfying 2Tr (A" A'') = 25"^ 

The second and third terms in ( |1.2D refer to the chiral quarks 

v|/R,i = {i(l±75)^^ f = l,---,Nj} 



are the SU{3)c 
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where Nf is the number of flavours, and each ip-^ is a color triplet. The quark mass matrix 
M, which is in general non-hermitian, is an Nf x Nf matrix in the flavour space. 
The second piece in ( |1.1| ) is given by 



Cg = - 9 A , where A = ^-^Tr 



:i.3) 



where < < 27r is an angle parameter, and the dual field strength tensor is defined by 
F =-e rF"^!^ 

For the coherence of the discussions it is convenient to specify the symmetries of the 
QCD Lagrangian first. The main difference between Co and Ce is that the latter does 
not respect the time-reversal (T=CP) invariance. This is clear from the fact that (in the 
temporal gauge: Aq = 0) the first term in Cq corresponds to — whereas Ce does to 
Ea ■ Ba, where Ea and Ba are the non-Abelian electric-like and magnetic-like fields. In 
fact, it is this CP-violating character of Cg which makes it fundamental for calculating 
the EDMs, which are inherently CP-odd quantities. 

Furthermore, since the gluons are flavour-diagonal, the massless quarks possess a 
UilNf) X Uji{Nf) global symmetry in the massless limit; M — * 0. In practice, compared 
to the dynamical scale of QCD {Aqcd ~ 300 MeV) one can take u, d and s quarks 
approximately massless. To this approximation the global symmetry of ( |1 . 1| ) becomes 
Uii^'i) X Ur{?i). The diagonal group ?7l+_r(3) ~ S'f/(3)y x U{1)b is also respected by 
the vacuum state since {i) hadrons obey an approximate 5'f/(3)y symmetry ( i.e. the 
Eightfold-Way)., and {ii) the baryon number is conserved. The remaining symmetries 
Ul~r{?>) ~ SU{?>)a X U{1)a are not manifest in the hadron degeneracies so that the 
dynamics should be such that the QCD vacuum breaks all these axial symmetries. Ac- 
cording to the Goldstone theorem, from the break-down of S'[/(3)a we expect eight pseudo 
scalars which can indeed be identified with the set [vr^, 7r°, vr^, , K^, K^, K~ ,r]]. How- 
ever, there is no pseudoscalar meson in the spectrum to be identified with the Goldstone 
boson of U{1)a- Hence, in the massless limit, there should be a remnant global symme- 
try of the QCD Lagrangian given by the phase rotation ip — >■ "^^ ipi- However, this 
symmetry is realized neither linearly (by the vacuum invariance) nor nonlinearly (by the 
existence of a light pseudoscalar). 

As will be detailed in the following sections, the resolution of the U{1)a problem and 
the appearance of the CP-violating Lagrangian Cg are both related to the nontrivial 
vacuum structure of QCD. Although there will be explicit quoting of the relevant work 



in the analyses below, still an exhaustive review of these problems can be found in ||T2 
together with [O . 



1.2 QCD Vacuum 

For a clear understanding of the QCD vacuum it is convenient to discuss the topolog- 
ical properties of the continuous functions. To study the topological properties of the 
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continuous functions one can divide them into homotopic classes; each class contains 
the functions that can be deformed continuously into each other. As a warm-up case 
one can consider the 5*^ mapping. Let A be the set of points {a} on a unit cir- 

cle S^, and let B be the set of unimodular complex numbers |e*^}. Then the function 
f-y (a) = e'**^"""'"'''^! realizes the mapping {a} C ^ {e*'^} C S^. For a fixed n, f-y (a) 
forms a homotopy class as 7 varies . More precisely, given /^^ and /^^ with 70 7^ 71 
then these two functions can be deformed continuously into each other by the homotopy 
F{a,x) = eH'^°'+(^-^ho+x'n)] -^j^ere x G {0,1}. One can visualize /7(a) as mapping n 
points of one circle into one point of the other circle in other words one winds around the 
latter by n times. Every homotopy class is characterized by a winding number 

f'^'^ da —i dL (a) ,^ 

n= / — ■'^^ ' (1.4) 

Jo 27r (a) da 

In particular, mapping with n = 1 is fy\a) = e*", and one can obtain the mapping for 



any winding number n by using yf\ 

Though the mapping above is useful for illustrating the homotopy classes, 

a physically relevant case occurs only when the symmetry group under concern has an 
SU{2) subgroup. Here the question arises: "How much it depends on the gauge group 
being SU{2) ?" The answer is: Firstly, if the gauge group is U{1), it is easy to see that 
every mapping of into U{1) is continuously deformable into the trivial mapping (all 
the S"^ is mapped into a single point). Thus for an abelian gauge theory there is no analog 



of winding number. Secondly, there is a theorem by Raoul Bott |T8| which states that 
any continuous mapping of 5''^ into G can be continuously deformed into a mapping into 
SU{2) subgroup of G, where G is a simple Lie Group. Thus everything which will be 
discovered for SU (2) will be true for an arbitrary simple Lie Group, in particular it will 
be true for SU{n). 

As will be clear below, the SU{2) group structure is needed to have a one-to-one cor- 
respondence with the spherical geometry at spatial infinity. Therefore, we now discuss 



the mapping - a mapping from the three-dimensional Euclidean space 

(parametrized by three Euler angles) into the SU{2) space (characterized by three pa- 
rameters) . In fact, we now seek for solutions to the classical SU{2) Yang-Mills theory in 
Euclidean space (|r|^ = xl + x ■ x) with finite action 



<Se = ^2 I d^x Tr 



;i.5) 



that is, the field strength tensor must vanish at infinity, 

F^,{x) — . (1.6) 

\r\~*oo 

Normally this condition is equivalent to vanishing of Afj_{x) at infinity. However, a van- 
ishing field configuration is equivalent to Ud^U~^ with U G SU{2) thanks to the gauge 
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transformation 



a; = UA^U-^ + Ud^U-^ . (1.7) 
Therefore, F^j_y{x) still vanishes at infinity ( |1.6| ) for in the pure gauge] 

A,{x) — u-'d,u. 

One notices that the points at infinity in the Euclidean space are three-spheres, S^. There- 
fore, the gauge transformation matrix U G SU (2) realizes a mapping from ^^(Euc. space) 
to the S'[/(2) (group space). More precisely, the manifold of the SU{2) group elements is 
topologically equivalent to the three-sphere S^. This can be seen from the fact that 
U G SU (2) can be expressed as 

U = e''^-" ^UQ + u-a , for |a| << 1. (1.8) 

Due to unitary character of the gauge transformations, WU = UW = 1 , it is clear that 
Uf^ spans a sphere 

ul + u ■ u = 1 

that is the group space of SU{2) is topologically equivalent to the unit sphere, S^. Fol- 
lowing Coleman [Q, the previous expression (|1.4| ) for the winding number can now be 
written as 



AuApAx 



d^x A (1.9) 



where A^ here are in the pure gauge discussed above, and the quantity A is defined in 
0| ). This is the expression for winding number corresponding to the SU{2) 



gauge theory. From this very expression of the winding number ( |1.9| ) it can be shown 
that the action ( |1.5[ ) is indeed finite, Se > ^nn^/g^, and the corresponding gauge field is 



either self-dual or antiself-dual, i.e., F^^y = ±F^i.. 

The finite-action solutions of the classical Euclidean Yang-Mills system ( |1.5| ) are 
instanton solutions ||21]] 

where 7 is an arbitrary scale parameter, and the corresponding gauge transformation 
matrix has the form ^ 

UAr) = — (xq + ix ■ a) . 

\r\ 

Thus, the gauge field is of finite extension, and as |f| — ^ cxd it approaches the pure gauge 
form, as expected. For this very solution the action integral has the value Sir'^/gg, that is, 
this solution corresponds to the homotopy class n = 1 hence the subscript "1" in Ui{x). 
Clearly the solution for the nth homotopy class can be obtained by compounding U hj n 
times: 

t/n(r) = (f/l(r-))^ 
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Suppose that there are two vacua belonging to different homotopy classes n and m. 
Then the vacuum to vacuum transition amplitude reads as 



(n|e-^^V) = J [dA^]n-me'' (1.10) 
In imaginary time (or Euclidean space), this becomes 

where Se = —iS is the Euclidean action. 

Therefore the intanton configuration corresponds to tunneling between the vacuum states 
having different winding numbers. The rough estimate above, based on the finiteness of 
the Eucleadean action, shows clearly that the effect here is inherently nonperturbative, 
that is, the transition rate is enhanced in the limit of large gauge coupling. Besides this, 
the rate is enhanced for neighboring vacua; \n — ml = 1. In the limit of small gauge 
coupling, the transition amplitude is diminished leaving system in any of the (equivalent) 
vacua with inherently perturbative dynamics. In summary we have infinite number of 
vacua each being characterized by its homotopy class or the winding number. Another 
point of primary importance is the behavior of these vacua under gauge transformations. 
The gauge transformation that changes to A^'^^ works as 

= U,d,U{' + U,A;U{' (1.12) 

so that the collection of the vacua {\n)} is not gauge-invariant at all: 

U-^\n) = |n + 1) . (1.13) 

Thus, these infinitely many vacua are connected by topologically nontrivial gauge trans- 
formations. 

Since the vacuum states |n) belonging to different homotopy classes (having different 
winding numbers) are separated by the energy barriers (as dictated by the tunneling 
amplitude in ( |1 . 1 1|) ) , it is clear that the true vacuum state will be a linear superposition 
of all these vacua. Similar to the situation with periodic potentials in quantum mechanics 
(whose true ground state is the Bloch wave), we define the true vacuum to be 

|e) = 5^e-^|n) (1.14) 
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which we name as the — vacuum" . As expected, in contrast to {|^)} the 6 — vacuum 
is gauge-invariant. This can be seen by computing 



—ind 



n) 



E 



n + l) 



;i.i5) 



Therefore, 9 labels the physically inequivalent sectors of the theory, and in each sector we 
can work physical processes in a gauge-invariant manner. The different 9 worlds do not 
communicate with each other unlike {|n)} which can be communicated via the quantum 
tunneling. Obviously, 9 is arbitrary, and thus, there is no a priori way of determining it 
from the theory. 

The amplitude for a transition between two 9 worlds in the presence of an external 
conserved current J is given by 



-iHt 



E 



m,n 



{m\e 



-iHt 



^-i{n-m)e ^im(e -0) 

m,n 

S(e' - 



E 



-iuO 



u=n—m 



n)j 



+JA) 



;i.i6) 



where in the last step we have introduced the effective action density 
L^ff = L + AL with AL = - 9u = - g^^ Cg 



(1.17) 



where Cg is already defined in ( |1.3|) . Stating in more explicit terms, the entire effect of 
the 9 vacuum can be taken into consideration by introducing the effective action 



d X Tt 



- 9 



47r 



d xTy 



:i.l8) 



where we have rescaled the gauge field as —>■ Qs Afj_, for convenience. It is clear that 
this effective action is equivalent to ( |1.2|) . 

In concluding this section we note that the QCD vacuum has a nontrivial structure 
described by the gauge field configuration corresponding to "pure gauge" . The finite ac- 
tion solutions of the Euclidean SU (2) Yang-Mills system (|1.2|) split into infinite number 
of equivalent subsets corresponding to the homotopy classes of S''^ (Euclidean space) — >• 
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S'^(SU(2) space) mapping. The gauge field configurations are instanton solutions in that 
they correspond to tunneling between the vacua with different winding numbers. There- 
fore, the Yang-Mills system is characterized by infinite number of vacua communicated by 
quantum tunneling. Moreover, the vacua are not gauge invariant; gauge transformation 
with winding number u transforms the vacuum state \n) to \n + u). 

The true ground state of the theory is given by a linear superposition of these vacua 
whereby defining the so-called 6 vacuum. The 6 vacuum is gauge-invariant, and each 
value of the parameter 9 corresponds to a physically distinct state of the system. This 
parameter is in general arbitrary the theory gives no way to determine it. 

There is no communication among distinct 6 vacua; moreover, the corresponding tran- 
sition amplitude is described by an effective action which differs from the original Yang- 



Mills action by in ( |1.3| ). The most interesting property of this additional piece is that 
it is odd under the CP transformation. Therefore, the appearance of Cg is due to the 
nontrivial nonperturbative topological structure of the QCD vacuum. 



1.3 Effective Vacuum Angle 



As was mentioned when discussing the symmetries of the QCD Lagrangian, there is 



an excess axial global symmetry {U{1)a) of ( |1.2|) which has no imprint in the hadron 
spectrum. A resolution of the U{1)a problem rests on the anomalous nature of this 
global symmetry and the nontrivial structure of the QCD vacuum discussed above. The 
conserved current, for massless u, d, and s quarks, reads as 



;i.i9) 



At tree approximation this current is conserved because its divergence, which is pro- 
portional to the quark masses, vanishes identically in the massless limit. However, this 
current is no longer conserved if one goes to one loop level 



(9 7^ 



2 Nf gl A 



;i.2o) 



where Nj = 3 for the present case. This result is independent of the looping quark masses. 
Moreover, this result is exact, that is, it does get no contribution from higher loops (the 
famous ABJ anomaly p2|, |2^). It is this quantum mechanical nonconservation of 



which forbids an exact U{1)a symmetry pj]. In this sense the QCD Lagrangian has no 
excess symmetry compared to the observed hadron spectrum, the partonic (quarks and 
gluons) and hadronic symmetry properties agree. 

Though this observation solves the U{1)a problem, there is more to be done with 
the QCD anomaly of J^. First, one observes that A can, in fact, be written as a total 
divergence 



levr^ A with K„ 



Tr 



F A 



'1.21) 
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with, however, the fact that is not gauge invariant at all. Despite this, one can embed 
the anomalous divergence in (|1.20|) to define the a new current 



'^5 



1^ 



2 Nf 



167r2 



;i.22) 



which is conserved in the massless quark limit. Therefore, the anomalous character of 
(that is the mismatch of the classical and quantum symmetries) is now removed via 
(|1.22|) . However, under a gauge transformation its conserved charge 



d^xJl 



;i.23) 



can be shown to shift as follows 



7° 

^5 



2NA—K') 



'Arc 



.f^-2Nfi^K') + 2N, 



f ' 



:i.24) 



where in arriving at the last line of the equation use has been made of ( |1.22| ). This is 
a very important result as it implies that the ^-vacuum is not invariant under chiral 
transformations. To see this one computes 



(1.25) 



which says that after a chiral rotation by an angle a the 6 parameter is shifted by 2Nfa. 
Therefore, similar to the fact that the gauge transformations shift the winding numbers, 
or equivalently, cause transitions between vacua of different homotopy classes, the chiral 
transformations shift the ^-vacuum depending on the number of massless quark flavors. 
This very shift of the 6 parameter under a chiral transformation has important implica- 
tions concerning the quark sector of ( [1.21 ). 

It has long been known that CP is not a respected symmetry in the neutral i^-meson 
(K^) system since the mixings as well as the decays of these mesons apparently violate 
the CP symmetry ?]. Since the QCD vacuum angle 6 does not contribute to CP, 
it is clear that there should be flnite CP violation in the electroweak sector (which will be 



throughly discussed in the next chapter), or simply, in the quark mass matrices in ( |1.2| ). 
The sensitivity of the ^-vacuum to the chiral transformations thus implies a hand-shaking 
of the vacuum angle 9 and the phases of the quark mass matrix M.g. 
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In general, the mass matrix of quarks is non-Hermitian, so that part of ( p..2|) containing 
the quark mass terms can be rewritten as 

Lm = i^R^MabiJL, + ^LSM^)abi^R, , (1.26) 

where a, b are flavour indices, and Mab are elements of Aiq. However, by a chiral SU{Nf) 
( here Nf = 3) transformation 

^RiL) = Ur^dMl). Ur^l) e SUr(^l){N), (1.27) 

one can diagonalize Aiq with positive elements up to an overall phase factor. Thus, 
without loss of generality the elements of A^, can be assumed to have the form 

Mab = mJab e'" , (1.28) 

where it is clear that the common phase factor here cannot be rotated away. Using this 
form of Aiq in ( p,.2| ) one observes that after a further chiral rotation, similar to ( |1.25| ), 



by the angle = —6/(2Nf) it is clear that the topological term Cq disappears from the 
QCD Lagrangian (|1 . 2|) . However, this very chiral transformation acts on the quarks as 
e~'^'^'''^°tp. Therefore, the expression for quark masses ( |1.28| ) now goes over to 



M,b = mJab e'" (1.29) 

where 

is the net phase of the quark mass matrix. One notices that now there is no 6 term in the 
original QCD Lagrangian ( |1.2| ); however, the same phase is now moved to the quark mass 
Lagrangian through the angle p. In general the hand-shaking between the CP phases in 
the QCD vacuum {6) and the phases in the quark mass matrices (p) is expressed through 
the effective 6 angle: 

e = ~Nfp = e - ATgDet[Mq] (1.31) 

where, obviously, Arg Det[A1g] = Nj p. 

Obviously it is the nature of ( |1.31|) which makes the CP violation by strong interactions 



nontrivial. Indeed, if 6 where stable under chiral transformations then there would not be 
any contribution to 6 from the phases of the quark masses. In such an instance, setting 
6 to some value (say, zero) would be a natural operation because 6 cannot anyhow be 
determined from the theory itself. However, even if one chooses 6 zero initially, this choice 
cannot be kept unchanging due to the nonvanishing (because of finite CP) phases of 
the quark mass matrices. Therefore, the CP violation in the QCD Lagrangian (|1.2|) is not 
removable at all. One can interpret it either as the topological Cq piece in ( p. . 2|) with real 
quark masses, or equivalently, as the quark mass Lagrangian with the phase ( |1.30 ) with 



£g absent. Except for this hand-shaking, there is an unremovable source of CP violation 
in the QCD Lagrangian with primarily nonperturbative, topological origin. In the next 
section we discuss an observable effect of 6: the EDM of the neutron. 
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1.4 Computing EDM of Neutron 

As was stated in the Introduction, the nonvanishing, unremovable 9 induces a finite EDM 
for the neutron. Calculating this effect with a high degree of precision is difficult as it is a 
long distance effect, and it may not be possible to employ the parton picture. In actuality, 
one must employ either an empirical model for the neutron structure (say, the Bag Model) 
[ p^ or make use of the current algebra techniques . However, the main purpose of the 
work here is to derive the effective Lagrangian {SCcp below) responsible for CP-violating 
effects in the strong interactions. Therefore, instead of a detailed calculation we will quote 
the results of the earlier works in estimating the neutron EDM. 

As a starting operation one should first determine the CP-violating Lagrangian which 
is responsible for developing the EDM. This Lagrangian, 6C{CV), follows from the quark 
mass Lagrangian. Consistent with all the above calculations we consider three light 
fiavours u, d and s. As suggested by ( p..31[ ) if 9 or at least one of the quark masses 
vanishes (that is, the determinant itself vanishes) then 6C{CP) — 0. To make it clear, 
consider an appropriate unitary transformation reducing the quark mass Lagrangian with 
the mass term in ( p..28| ) to 



Cj^-j — > C\-^ = {e"-'')muURUL + rUddRdi + rrisSRSL + H.C. (1.32) 

where we have arbitrarily transferred the entire phase content to the u quark mass. From 
this we can see that if m^i then p drops out. Therefore, 5C{CV) cannot be simply part 
of the quark mass Lagrangian having imaginary parts. As mentioned in the introductory 
section, the spontaneous break-down of the axial SU{'i)A symmetry of the QCD La- 
grangian produces a total of eight Goldstone bosons. In this case the vacuum is infinitely 
degenerate, and the naive perturbation theory is not reliable. To determine 5C{CV) we 
have to select one vacuum state out of this infinitely many sets, which of course leaves 
the physical content of the theory unchanged. To leading approximation the vacuum is 
fiavour-blind. The correct perturbation 5C should be chosen unitary-equivalent quark 
mass matrices ( p..27| ). Moreover, this perturbation should cause a minimal shift of the 
vacuum energy with no danger of destabilizing it: 

iSimn)= min {n\Cj^,{UR,Lm. (1-33) 

where Vl stands for the vacuum. This constrained selection of states with unitarily equiv- 
alent quark mass matrices can be done with the help of Dashen's theorem ([Q which 
assumes that {i) the transformations should be pure chiral, = Ul, and (ii) 6L is 
fiavour-blind. With these conditions at hand, one can parameterize 6L as 

6C = i^RiiJa + iu)%l)le-'^ + H.C, (1.34) 

where and uo are assumed to be real parameters and are fixed by 
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2 2,2 



(1.35) 
(1.36) 

when /id > 0. (1-37) 
Now we will derive the correct form oi CP violating perturbation. Consider ( |1.34|) 

= [jJa{L + R) + iuj{L - R)] tfj^e'^ 



where we have used, 



L 



1+75 



R 



1-75 



2 ' "2 
Hence 6C{CP) is the CP nonconserving part we need excluding 6. Therefore, 

SC{CP) = T0J^il5'^ 



where uj can be solved by inverting (|1.35| ) as: 



e = -Are 
3 ^ 



(/ii + iuj){jj2 + ioo){iJ.3 + iuj) 



:i.38) 



;i.39) 



;i.4o) 



;i.4i) 



therefore 



tan^ 



-tan Arg 
3 ^ 



U) \ UJ \ LjJ 



^ /^l ~^ ' ' M1A^2M3 

3 l + u^H^ + ^+ 1 



^ ^l\^i2 Mi/is • 



;i.42) 



In the limit \6\ « 1, and when also uj is small ( in order to neglect u; , a; ) we obtain 

/^lyU2/i3 



= 3 



yUi/i2 + yUl/i3 + /^2/i3 



-d 



1.43) 



From ( |1.35| ) for small uj we can write = yU^, which upon substitution in the last 
expression yields, 



= 3 



171117127713 



mim2 + mim^ + m2in3 



-e 



;i.44) 
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which can be used in (|1.4CI| ) to rewrite it as 

S/:{CP) = T th^a^ ^(?Z75^). (1.45) 

m„mrf + munis + nidms 

Derivation of this Lagrangian is the most important step in studying the CP-violating 
effects in the strong interactions. It is this very Lagrangian with which there arise CP- 
violating effects in the hadronic interactions. The EDM of the neutron is one such fun- 
damental quantity through which the the parameter becomes observable. Below we 
will sketch its derivation without going into details of the hadronic model adopted. Two 



existing calculations due to |jT^ and ||25| are in good agreement. 

We now choose to sketch the derivation in [^. The neutron EDM is defined via the 
correlator 

in + ^ (k) I j d'^x 6CiCP)\n (^p - ) = - D„ (A;^) ua^^k'u , (1.46) 

where J^™' electromagnetic current of the quarks. The value of the formfactor D{k'^) 
for on-shell photon (that is, D(/c^ = 0)) is the EDM of the neutron. The Figure below 
illustrates the correlator above. Here the black blobs designate the usual n'np and 7r^7r+7 
vertices. The grey blob, however, shows the CP-violating nnp vertex generated by SCcp 
in the correlator above. 




Figure 1.1: The mechanism of generating the neutron EDM. The black blobs are the usual 
CP-conserving hadronic couplings. The grey blob shows the CP-violating n^np vertex 
generated by 5Ccp- 



In general the calculation of the correlator ( |1.46| ) involves a summation over all inter- 
mediate states 

Y,{n\rr\^) X {X\5C{CP)\n) 

X 

where |X) = {|A^), |vrA^), IvrvrA^), ■ ■ ■ }. However, as shown in [^], it is the IvrA^) contri- 
bution that dominates the amplitude, in accordance with the figure above. The reason for 
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the dominance of this particular intermediate state follows from the fact that the neutron 
disassociates into charged constituents (needed to emit a photon) in acquiring an EDM. 
The largest contribution comes from the lightest constituents ( pion) because it attains 
the greatest distance from the center of the charge distribution. 

Due to the explicit CP violation in the system the pion-nucleon interactions are now 
generalized to 

^nNN = ■ Na {ig^NNl5 + W^) N (1.47) 
where gnNN reflects the CP violation. A direct evaluation of {7i°'N\6Ccp\N) leads one to 



g-KNN 



g^NN = - Ox X ^ ■ 1.48 

which is proportional to 6 as expected. The remaining vertex in the correlator ( p,.46| ) is 
J^'^\ttN) is known to be proportional to 'y^cr^^k'^ [ p^ . 
Then a direct computation of the diagram in the figure above gives 

rf..D„(P = 0) = ^^^^lnf^') , (1.49) 



which equals numerically ~ ^ x 2 x lO^^^e — cm ^5 . 

Here one notices that this numerical estimate is quite close to the results reported by 
Baluni WM though two calculations adopt different methods in computing the hadronic 



matrix elements. However, a more interesting coincidence comes by the recent estimate 
29l of dn using the QCD sum rules [^. In this calculation the main object is the two- 



point correlator of the neutron current (that current which excites neutron from the QCD 
vacuum) in that background containing the electromagnetic field and the ^-term derived 
above. 



1.5 Remarks 



The numerical value of the neutron EDM given above is unfortunately far above the 
experimental upper bound: \dn\ ~ lO^^^e — cm ||2^. This means that the effective QCD 
vacuum angle should satisfy \9\ < 10"^ ! This requires a huge fine tuning of the pure 
QCD angle 6 and the phases of the quark mass matrices (|1.31| ). That one has such a 



small value for 6 instead of the expected order of unity poses the well-known strong CP 



problem - a CP hierarchy, or a fine-tuning, or a naturalness problem |31 



We will not attempt to discuss possible solutions to this issue instead we will summarize 
shortly existing proposals. There are two main ideas towards a solution to the strong CP 
problem: Relaxation and cancellation. The former refers to the celebrated Peccei-Quinn 
solution (see also reviews |33|) according to which 6 (being now a dynamical variable) 
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relaxes to zero and remains so to all orders in perturbation theory. More explicitly, one one 
(i) promotes the phases of the quark mass matrices to dynamical variables (corresponding 
to the Goldstone bosons of spontaneously broken global symmetries), then {ii) computes 
the instanton-induced effective potential for 9, and finally (Hi) shows that this potential 
is actually minimized for = 0. Therefore, in this picture, 6 relaxes to zero dynamically 
resulting in a vanishing EDM for the neutron. This Peccei-Quinn idea has generated 



several other versions in quest for an experimentally viable model: Weinberg and 
Wilczek put the Peccei-Quinn scale to the weak scale which disagreed with the 



phenomenological constraints then, later Dine et.al shifted it to the unification scale [36 



As another version, instead of working with the quark mass matrices directly, Kim and 
Shifman et.al introduced additional heavy color triplets [^. In all versions of the Peccei- 
Quinn idea the effective vacuum angle 6 possesses a dynamical character and then relaxes 
to a purely CP-conserving point via the background instanton effects. 

The other proposal to solve the strong CP problem runs via an appropriate choice of 
the quark mass matrix Aiq (which is block diagonal in up- and down-sector quark mass 
matrices ) such that Det[A1„A1d] is real whereas the CKM matrix (discussed in the next 
chapter) keeps having a finite phase , as originally proposed by Ann Nelson. With such 
a cancellation of the phases coming from the quark mass matrices, the interplay between 
the original QCD vacuum angle and that of the quark sector is lost, so that one can take 
6 = from the scratch. For instance, in the Babu-Mohapatra model ^9\, the theory is 
left-right symmetric (parity-conserving) at higher energies so that the topological term 
is not allowed to contribute at all. Subsequently a finite 9 could be generated were not it 
for the flavour structure of the quark mass matrices. 

At this point it is convenient to comment on the situation in the supersymmetric 
models concerning both ideas above. The supersymmetric models provide enough global 
symmetries replacing the Peccei-Quinn symmetry together with novel sources of CP vi- 
olation coming from the soft terms [0. The Peccei-Quinn solution can therefore be 
adapted to supersymmetry for not only solving the strong CP problem but also solving 
its own hierarchy problems p] , Finally, the Peccei-Quinn idea can be generalized 
also for examining the CP hierarchy problem of supersymmetry in a dynamical way 
Apart from the operation of the Peccei-Quinn mechanism in supersymmetry, it is possible 
to implement supersymmetric falvour models which incorporate the Nelson-Barr models 
^ which can be particularly useful for having observable supersymmetric CP violation 

ni- 



ls 



Chapter 2 

The Neutron EDM: Electroweak 
Interactions 



In the framework of the Standard Model (SM) of electroweak and strong interactions 
the other source of CP violation is the single complex phase in the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix [|l^, |16]. To examine it's consequences for EDM of the neutron 



we will first briefly review the electroweak sector of the SM with particular emphasis on 
the CP-violation currents and parameterizations of the CKM matrix. 

2.1 Electroweak Lagrangian and the CKM matrix 

The particle spectrum of the SM consists of gauge and Higgs bosons together with three 

families of quarks and leptons (See [?] for a review): 

Fermions 

Leptons: (^""^ ^ en,m,r^ 

Quarks : ( ^ ^ , f ! 1 , ( [ ] ur, Jr, cr, 
Gauge Bosons 



d ' \ S ' V ^ / T 



Scalars 



( iv; w'^ wi ) i B„ 



*= !o . !.o (2.1) 



The interaction between the fermions and the gauge bosons has the form: 

3 6 

£(/, W,B) = Y, + I'bP^'r + Wl] + E ^BP'<tB. (2.2) 

j=l i=l 
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where j is the family index, 

" ( ) ' = ( ) ' ^R = '^R^ ^R = dR,--- 



L \ / L 

and the two covariant derivatives are defined as: 



D'^ = d,~igr^B,. (2.3) 

where gi and g2 are respectively, the f/(l)y and SU{2), coupling constants, and (T'*(a = 
1, 2, 3) are the (weak) isospin matrices. The interaction term £(/, (p) involving the fermions 
and the Higgs fields has the Yukawa form, 

3 3 

^(/>'^) = E{(^^)/«} + E {(^;).#4+(/^.wi<^'^4| (2.4) 
j=i j,k=i ^ 

where the charge-conjugated scalar doublet reads as 

= i(T2(, 

with both $ and $^ transforming as a (weak) isospin doublet with opposite hypercharges. 
After the spontaneous symmetry breaking SU{2) (8> U{1) — > U{1)em, the gauge bosons, 
fermions and the neutral scalar field, , acquire non-zero masses through the Higgs 
mechanism. After develops a vacuum expectation value v, one can expand the Higgs 
doublet around this ground state so that (p^) becomes: 



j=i ^ ^ ^ j,k=i ^ ^ ^ ' 

(2.5) 

where the fermions mass matrices (in the fiavour space) are defined by 



V 



{mjk)D = -{h'g)jk^ , (2.6) 

with {mjk)u and {rnjk)D being 3x3 mass matrices for up and down quarks, respectively. 
In order to write the Lagrangian in terms of quark mass eigenstates, the mass matrices 
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{fTijk)u and {rrtjk) d have to be diagonahzed. This can be done with the help of two unitary 
matrices usually denoted by V^^ and V^'^ (similarly for down quarks): 

Vl^muV]^^'^ = {mdiag)u = Diag.{mu,mc,mt) 

"-"mD^^""""^ = {md^ag.)D = Diag.{md, m„ m,) (2.7) 
with V^^'^v^ = 1 , etc. Considering only up quarks we can write the mass term as: 

= ^I^{md.ag.)u{V7^R). (2.8) 
which shows that the physical quark states are: 



^Phys ^ yup^^ ^ yup 




^Phys ^ ydown^^ ^ ydown ^2.9) 



Then in terms of the mass eigenstates, ( p.5| ) can be rewritten as: 

6 3 



The identification of the parameters mi- , rriq. with the lepton and quark masses is now 
clear. Since we have written the term £(/, 0) in terms of the physical quark fields, we will 
express other terms of the Electroweak Lagrangian in terms of the physical quark fields 
like£(/,iy,5). 

Considering first the 'neutral current', we see that neutral current part of C{f, W, B) is 
manifestly flavor diagonal. Written in terms of physical boson (VF^, A^) and fermion 
fields: 



2V2sin2 e 



-Z^{hL - Q sin^ ew)i + eA^Qi 



w 



U (2.11) 



The neutral current interaction induced by the Z-exchange violates P and C but conserves 
CP. 

It is important to emphasize here that the Higgs-fermion Yukawa couplings are flavour 
diagonal so there are no flavour changing neutral currents (FCNC). Thus, all the flavour 
changing transitions in the Standard Model are confined to the charged current (CC) 
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sector. From ( |2.2| ), concentrating only on the quark sector, we obtain the charged current 
interaction lagrangian, 

3 

C^^ = ^ y y'l^^W+di + h.c. (2.12) 

2v^sin2^-^ L/ L 



w 



i=l 



Making use of ( pl9| ) in ( p.ll| ) we finally get: 



^cc ^ % Yl ^i^rV2'viy;^/°""Vf + h.c 

2V2sin .^^ 

= E i<'n'rw;iivrvt'"''%id^,''y + h.c. (2.13) 

2V2 sm Ow 

Therefore, the charged current J^*" which couples to the is, 

J^^^ = ( c, t) ^,VcKM s j (2.14) 

rup-rrdown] 



where Vckm = Vl^'Vl"'"''^ is the CKM matrix 0, 0. It is a generalization of the 



Cabibbo rotation [|TB| to three quark flavours, and it was introduced to keep the unitarity 



of weak interactions. The charged current Lagrangian has a (V-A) structure, hence it 
violates P and C maximally. 

In general, CP^ violates CP due to the possibility of some non-trivial phase in Vckm- For 
instance, consider (1.12), dropping the superscript Phys and using the CP properties of 
fermionic flelds we get: 

3 

2V2 sm 

^ 47'^iy-(W/).,ul + ul7'^l^;(V^ci.M)..^i (2.15) 

where we have used W'!^ — > ~W~ . 

^ fj, 

Thus the CP conservation requires the matrix V to be real. It is real for two fami- 
lies, and hence, for two families of quarks CP is automatically conserved. But in case of 
three families of quarks, as was suggested by Kobayashi and Maskawa, the CP is violated 
due to the presence of complex phase in Vckm- 

Vckm being a unitary matrix satisfies, V^j^j^^^Vckm = 1- The matrix elements of Vckm 
are determined by the charged current coupling to the bosons. Symbolically this 
matrix can be written as: 

Vud Vus Vub 

Vckm = \ Kfe | (2.16) 

Vd Vs Vb 
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As long as the unitarity is maintained, one can parameterize the CKM matrix using 
three Euler angles (appropriate for 3 x 3 matrices) and a complex phase Sqkm- There 
are different parameterizations of CKM matrix. The original one, due to Kobayashi and 



Maskawa [T^ was constructed from the rotation matrices in the flavour space involving 
the angles 9i{i = 1, 2, 3) and the phase 5c km, 

VkM = i?23(^3, 5cKM)Rl2{ei, 0)i?23(^2, 0) (2.17) 

where < 6*^ < |, < Sckm < 27r, and Rij{0, 0) denotes a unitary rotation in the 
plane by the angle 9 and the phase (p. Then a possible representation for Vkm is: 

/ Ci -S1C3 -S1S3 \ 

Vkm = C1C2C3 - S2Sse'^o^^^ C1C2S3 + 32036^^^7*-./ (2.I8) 

\ S1S2 C1S2C3 + €2836'^^^^^ C1S2S3 - C2C3e^^C^^" / 

with Ci = cos 9i, Si = sin 9i. This reduces to the usual Cabibo form for 6*2 = ^^3 = with 
9i9c identified with the Cabibo angle. 

However, in discussing several flavour-changing processes it is useful to consider an ap- 
proximate form of Vkm- Indeed, Wolfenstein |^6|, has made made the important obser- 



vation that (empirically) the Vkm can be expressed as 

\Vii\c^ 1, 1 = 1... 3 
\Vi2\ ^ IV21I ~ A 
IV23I ^ IV32I ~ A^ 

IV13I ~ \V3i\ ~ A^ (2.19) 

with A = sin^c = 0.221. With this hierarchy of elements in terms of the Cabibo angle 
one can write 

/ l-iA2 A AX'{p-ir])\ 

Vwolfenstein = ( "A 1 — |A^ AX'^ j (2.20) 

\AX%l-p-ir]) -AA2 1 / 

where A ~ p ~ 1 and \ri\ << 1. This form of the CKM matrix is rather common so that 
the improvements in the experimental determinations are generally expressed in terms 
of the parameters here. Having briefly reviewed the particle content of the Electroweak 
Theory and the CKM matrix we are now in a position to describe the consequences of 
CP violation in Electro-weak sector for the EDM of neutron. 



2.2 Computation of Neutron EDM and Form Factors 

In the sections which follow we will describe the calculation of the neutron EDM in the 
electroweak sector. It is clear that the end result will be proportional to Sqkm, the only 
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source of CP violation in the electroweak theory. In estimating the neutron EDM we 
adopt the quark model where the neutron is composed of three valence quarks u, dm and 
d. In principle, if one calculates the EDM of a free quark q at the W-boson mass level 
then it is a straightforward issue to move it down to the nucleon mass level using the 
appropriate renormalization group equations (RGE). Then a rough expression for dn in 
terms of the quark EDM's reads as 



dn ~ \ {dd)MM - I (du)MM (2-21) 



where the subscript Mn shows that these moments are calculated at the nucleon mass 
scale. In estimating the contributions of the physics beyond the SM [e.g. supersymmetry 
48] or two-scalar doublet models this relation works well, and puts stringent bounds 



on the CP violation sources in the underlying model. Here one notices that the unknown 
parameters of the two-doublet models and the CP hierarchy problem in the supersymme- 
try can be avoided (at least partially) in models where supersymmetry is broken around a 
TeV leaving a two-scalar doublet model with less unknowns and appropriate CP violation 
sources 



50, 40 



u q 



(a) (b) 



Figure 2.1: Diagrams for induced u-quark dipole moment :the photon line is to be attached 
to all the charged lines in (a) one-loop and (b) two-loop diagrams. 



However, in estimating the contributions of Sckm to the neutron EDM one faces 
with certain difficulties associated with the EDM of a single free quark. After a detailed 
repetition of Ellis-Gaillard-Nanopoulos ^T|] analysis Shabalin |Q has shown that there 



is no contribution to a quark EDM up to three loops. This arises mainly from the 
cancellations among different diagrams due to the unitarity of the CKM matrix. To 
illustrate the situation it is convenient to refer to Fig. 2. 1. The one-loop graph (graph 
(a) here) in this figure is a self-conjugated graph which cannot contribute to CP violation. 
At the two loop level, there is a possibility of obtaining "non-self" conjugate diagrams, 
as seen from Fig. 2. 1. (b). Indeed, this diagram yields a vertex of the form 

edn{k'^)uaf,^k''-f5uAf' (2.22) 
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where dn{k'^) in the hmit — > is the EDM, and is suppressed by the usual GIM factor 
of Am^/M^ = ^VqiV*,.m'^M^^. Unhke the expectations of it was shown by |^ 



that these two-loop contributions sum up to zero. Hence one is left with the three-loop 
and higher contributions for the EDM of a single free quark though even the three-loop 
contributions were already shown to partially cancel ||53|. These cancellations bring the 
d-quark contribution down to 10~^^ e-cm. 

After illustrating the smallness of the single quark EDM's with electroweak theory, we 
note that the problem here can be circumvented through a more realistic approach to the 
problem. Indeed, one has to take into account the multi-parton content of the neutron 
in calculating the EDM, that is, instead of dealing with a single free quark, one should 
consider all three quarks simultaneously. The original analysis by Nanopoulos-Yildiz-Cox 
17 1 has thus dealt mainly with the diquark moments, where the third valence quark was 



taken as the "spectator" of the process in charge of balancing the spin, charge and the 
kinematics of the nucleon. Therefore, below we shall be analyzing the neutron EDM (in 
parton language) taking into account two (very weakly interacting) quarks into account. 
It will be seen that with such a multi-parton language the nucleon's moment turns out 
to be essentially a one-loop effect. 

2.2.1 Diquark Interaction and Neutron Moment 




(a) 

Figure 2.2: Type of diagram contributing to the EDM. Quark lines are labelled according 
to leading order contribution; external lines are labelled with momenta {Q << P, P', k). 



Depicted in Fig. 2. 2 are the class of diagrams contributing to the EDM of a two- 
quark system where the interaction between the two quarks is "weak", that is, mediated 
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by the boson. The shaded blob in this figure is the quark-quark-photon vertex which 
can get contributions from three possibihties shown in Fig. 2. 3. One keeps in mind that 
and 7 can be attached to any of the four external legs of the quarks. Schematically, 



y 
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w 
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(c) 



Figure 2.3: The three contributions to the three-particle vertex of p. 2.1 



the amplitudes in Fig. 2. 3 add to give a transition amplitude of the form 

(a form factor) x -^1^0) Pl 

where is the quark momentum, and is the momentum of photon of which the 
quadratic terms have been neglected. Direct compoutation yelds a five-particle amplitude 

r^,, = {{SGF/V2)Mlf V F{mlml)d{P-k) 
X P,^Z^^l±^,^P,u[P'-k) 

X A,A(P-P'MP' + %Virf(P + fc), (2.23) 

where A^a is the propagator in the unitary gauge, and V = (y'^)dtVts(y'^)suVud is 
a reparametrization invariant combination of the relevant elements of the CKM matrix. 
The quantity F is the form factor, which will be detailed below. 

The EDM is given by the coefficient of F^^^ ~ QA'^ — Q^A, so only T^^^ needs be deter- 
mined for calculating the diquark EDM. To continue the analysis we need to determine 
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first the quark spinors in the diquark rest frame. Parameterizing 7^ as 

(2.24) 



1 \ / a 
-1 j ' I -ff 



the quark spinor q = (u, or, d) can be written (as two-component spin states in non- 
relativistic approximation) 

q(p) = [2mip^ + m)]-'/' (^P' ^J^ ^ (2.25) 
where Qa = \ ] or (?) are the usual two-component basis vectors. Using 



X 7 V 1 . 

this construction in the rest frame of the diquark together with an averaging over the 
orientation of relative momentum of the quarks, we obtain the mean value of the five- 
particle vertex ( p.23|) as 

In deriving this equation we have made use of the spherical symmetry of the nucleon 
wavefunction so as to have the neutron spin function 

1^, ±) = ±^ {l±)u (I + -) + I - +))dd - 2|T)u| ± ±)dd} (2.27) 

after using the appropriate Clebsch-Gordon coefficients. Between such states the evalua- 
tion of the spin-dependent operators can be done using 

Q--l-«-l = |J-?"'f'" (2.28) 
I +Q ■ , for d 

where <?„ are the Pauli matrices for the neutron spin, defined in the basis ( 2.27 ). Another 
quantity in ( p.26| ) to be mentioned is = Efp' / m involves energy, momentum and mass 
of either quark in the diquark rest frame (remember we are using m„ = md = m). 

In addition to Fig. 2. 2. there is another diagram with the same structure, however, 
occuring at the other d quark leg. This diagram gives the same result except for the fact 
that V is now replaced by the V*, and the overall sign is changed. So the net result will 
be given by ( p.26|) with V replaced hj 2 i ImV, i.e. 

(PL) 4 (^MSy 2n,.V ^ ±-KQ.,. . (2.29) 

— * 

Integrating this result over the quark phase space, supplying a factor of 2 for the d quarks, 
we determine the appropriate coefficient to be identified with the EDM of the neutron 
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where A = {mq\Pq\^) in the neutron state to be evaluated in the diquark rest frame. 
One notes that A ^ for both iriq (ultrareletivistic hmit) and m„ — Sm^ 
(nonrelativistic hmit). Assuming that averaging in the neutron state is given by the 
phase space alone, A is a bounded function of irtq whose maximum being in between the 
two asymptotics lets one to make the rough estimate A < 10^'^ m^. Now we discuss the 
behaviour of the form factor F . It was introduced in (|2.23|) , and reads as 

F(m^^,m;) = - — — TT ie{ ^ e ^ + -e 2- * 



(47r)2 2M^ [ \2 M^J 2 V J 3 V 



where the three terms in the bracket correspond to the diagrams of Fig. 2. 3. From 
the experimental values of the top quark and mass we see that the terms in the 
parenthesis equals numerically ~ 10. It is clear that, if the top quark were not heavy the 
ratio m^/My^ would have brought about a GIM suppression factor, as usual. 

Using the parameterization of the CKM matrix introduced in the previous section we 
can write 

ImV = C1S1C2S2C3S3 sin 6ckm, (2.32) 

with Ci = cos 9i, Si = sin 9i, and the angles 9i and 5c km are the parameters of CKM 
matrix. 

Substituting ( p.31j ) into ( |2.30| ), the neutron EDM turs out to have the following explicit 
expression 

2 /8Gf\ ,„., ,„ -i 19 1 Mw m„ . 

The reparametrization invariant combinations of the CKM matrix elements are also 
contained by the amount of CP violation in mixings and decays of the and systems 



53 . For instance, assuming no new physics contribution, the experimental bound on the 



— mixing parameter requires |]TT] 



Cislc2S2C^s-iSin5cKM ~ 10 ^ (2.34) 

Using this parameterization, and supplying the masses of the quarks, we arrive at a 
numerical estimate for the neutron EDM : 

dn ~ 10-32 e - cm. (2.35) 

It is clear that this number is seven orders of magnitude below the experimental upper 
limit. It is plausible that the existing bounds will not change drastically unless there 
arises completely new high-precision experiments. Therefore, this gap of seven orders of 
magnitude is probably going to be closed by the physics beyond the electroweak theory. 
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2.3 Remarks 



In this chapter we have summarized the calculation of the neutron EDM when only the 
phase of the CKM matrix is present. The calculation involves several rough estimates 
for the hadron wavefunctions so there is always a certain amount of uncertainty in the 
theoretical prediction. However, the numerical result is seven orders of magnitude smaller 
than the present experimental upper bound, and this gap is unlikely to be closed by the 
theoretical uncertainties. Even if one (accidentally) forgets about the kaon system CP 
violation strength ( |2.34| ) the result is still 2-3 orders of magnitude below the experimental 
upper bound. It is in this sense that the EDM in the electroweak theory is much smaller 
than the (long-living) upper bound. 

Apart from two-doublet models |^) and supersymmetry [Q, in close similarity to 
what has been done in this section, one can discuss the neutron EDM in extensions of 
the SM with more than three generations, say four generations. In this case there are 
three phases in the corresponding CKM matrix with more reparametrization invariant 
combinations. However, there is still strong constraints from and systems ||5^ so 
that it is unlikely that the neutron EDM will be levelled to the experimental upper bound 
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Conclusion 



In this work we have presented a brief review of the neutron EDM calculation in the SM. In 
doing this we have analyzed the contributions of the strong and electroweak interactions 
separately. It is clear that the two force laws imply diversely different values for the 
neutron EDM: 

(dnjstrong I (dn) exp ~ > and (dn) electroweak I (dn) exp ~ ^0"^ • 

Each interaction here has its own source for violating the CP symmetry. These sources 
have nothing in common concerning their nature and function in the theory. That of 
the strong interactions arises from the nontriviality of the QCD vacuum concerning its 
topological and nonperturbative structure. However, the CP violation source of the elec- 
troweak theory is related to the number of families (for two families no source of CP 
violation, for four families there are two more sources), and the strength and hierarchy of 
the intergenerational mixings. 

The importance of these two contributions lies in the fact that they will always be 
present irrespective of what kind of extension of the SM is considered. For instance, in 



the supersymmetric models there are new sources of CP violation PD[ which generally 



exceed the experimental bounds by three orders of magnitude |HSI]. However, it is still 



meaningless to discuss the supersymmetric contribution alone, as the QCD contribution 
is there to exceed the bounds by several orders of magnitude unless a Peccei-Quinn type 
scenario is exploited (such as ETI 



or 

The EDM of the neutron receives contributions from strong as well as electroweak 
interactions as mentioned above. This remains true in any extension of the SM model. 
However, being devoid of any color degrees of freedom, the electron EDM receives no 
contribution from the QCD angle. Indeed, it is known that the experimental bounds on the 
neutron and electron EDM's differ by an order of magnitude as could be taken into account 
by their mass ratios; m^/me ~ 10. This expectation is confirmed by the predictions of 



supersymmetry |^8|] and two-doublet models It is here that one observes the "excess' 



nature of the QCD contribution which spoils the existing bounds by several orders of 
magnitude. For any underlying model, consistency among the EDM's and CP violation 
in meson systems is prerequisite for any meaningful prediction. 
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